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Abstract 

We investigate the relativistic kinematics of the spin-polarized colhsion of deuterium incident 
on tritium, producing ^He and a neutron. Within the context of special relativity, we apply the 
conservation of four momentum and the conservation of intrinsic spin, which leads to a system of 
ten equations. We impose initial conditions such that the deuterium is moving along the x-axis, the 
tritium is stationary at the origin of coordinates, and the classical spin vector of the deuterium (spin 
magnitude = 1) is along the +z-axis, while the classical spin vector for tritium (spin magnitude = 
1/2) is along the — z-axis. We expand the ten conservation equations to second order in velocities 
and we solve them for the velocity components of the ncTitron, its unit-spin-orientation vector, and 
the velocity components of the ^He nucleus, as a function of the incident deuterium energy. We 
find that this analytic solution agrees closely with the numerical solution of the ten (unexpanded) 
equations. For a given energy of deuterium, we find that there are two solutions, each solution 
having a unique velocity for the emitted neutron and helium nucleus. The two solutions are related 
to each other by reflection in the plane perpendicular to the deuterium spin and containing the 
initial deuterium velocity vector. 
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I. INTRODUCTION 



Sources of neutrons are currently used to excite nuclei in an unknown sample of material 
in order to determine its elemental composition. The technique is called neutron activation 
analysis, where neutrons are incident on target nuclei and a unique spectrum of gamma 
rays are emitted for each element that is present in the sample [1]. In this technique, mea- 
surements are made of emitted gamma rays that are either emitted almost instantaneously 
(prompt gamma-ray), or, that are delayed. The prompt gamma-rays that are emitted almost 
instantaneously come from a compound nucleus that is formed when the neutron is captured 
by the target nucleus in the sample. The technique is called prompt gamma-ray neutron 
activation analysis (PGNAA). When delayed gamma rays are measured, they are emitted 
by the decay of a radioactive intermediate state formed by the neutron irradiation. In this 
case, the technique is called delayed gamma-ray neutron activation analysis (DGNAA). In 
both cases, analysis of the composite emitted gamma ray spectrum often allows a precise 
determination of the elemental content of the sample [1] . A critical component in this type 
of detection is the source of neutrons. Ideally, we would want a narrow beam of neutrons 
that can be directed at the sample so that the damaging effects of neutron radiation are 
limited to a small solid angle. A reaction that is commonly used for neutron production is 
the collision of deuterium with tritium [2], 

d + t ^n + ^He (1) 

However, with unpolarized deuterium and tritium this reaction generally produces an 
isotropic distribution of emitted neutrons, as opposed to the narrow beam mentioned above. 
If the deuterium and tritium nuclei were spin polarized before collision, then the symmetry 
of the coUision, T {d , r?) '^He, would be reduced, and we would expect that there may be 
preferential directions for the emitted neutrons. Quantum mechanical calculations of this 
collision process are quite comphcated [2, 4]. However, it is possible to use a considerably 
simpler approach to analyze some aspects of the kinematics of the spin-polarized deuterium 
and spin-polarized tritium collision. We use conservation of four-vector momentum and 
conservation of intrinsic spin tensor, within the context of special relativity, to investigate 
the kinematics of this collision. 

This reaction in Eq.(l) has the advantage that it has a large cross section, which peaks 
at 5 b [2, 3] for a deuterium energy of 107 keV. It is beheved that at low-energy the reaction, 
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T(d,n)^i7e, proceeds via a short-lived intermediate resonant state of ^He with a spin angular 
momentum J = 3/2, and that this state quickly decays into a neutron and a ^He nucleus. 
At low energies, the reaction proceeds via zero orbital angular momentum, i = 0, so that 
the total spin J = 3/2 must arise from intrinsic spin of the rcactants [4]. Therefore, at low 
energies, the deuterium nuclear spin = 1, and tritium nuclear spin, St = 1/2, must be 
aligned to produce a total angular momentum J = 3/2 in the initial state. The products 
of the reaction are a neutron, with nuclear spin s„ = 1/2, and '^He which has nuclear spin 
Sh — 0. The total intrinsic spin of the product state (helium and neutron) is s = 1/2, so the 
rest of the angular momentum, the quantity J = 3/2 — 1/2 = 1, must be carried off as orbital 
angular momentum of the product state. Alternatively, there is also some small amplitude 
for the reaction to proceed, with zero orbital angular momentum, when the deuterium and 
tritium spins are anti-aligned. In this case, there is zero orbital angular momentum carried 
off by the products. We consider this latter case in the work below. 

The formalism we use was initially applied by J. L. Synge to discuss elastic collisions 
(without transmutation of mass) of particles having intrinsic spin [5]. Here, we use Synge's 
formalism to analyze the fundamentally inelastic collision, T [d , 'n) *He, when the reac- 
tants have a definite state of spin, in the sense of special relativity, where spin is represented 
as a four-tensor. The treatment we give below assumes that the deuterium has a zero impact 
parameter on the tritium, and therefore this limits our analysis to the case where the spin 
of deuterium and tritium are anti-aligned, i.e., the spins are pointing in opposite directions. 
Unfortunately, the condition of zero impact parameter is not easily relaxed because of com- 
plications in special relativity associated with the speed of propagation of interactions and 
the requirement of relativistic invariance [5]. As described above, at low energies, the zero 
impact parameter (zero orbital angular momentum in initial state) is probably the most 
important in the collision process since the reaction T {d , li) "^He will only occur when the 
deuterium and tritium are in close spatial proximity, due to the short range nature of the 
nuclear forces. 

Therefore, our goal here is to analyze the dependence of the velocity components (angular 
dependence) of the emitted neutron and helium nuclei, on the initial conditions, which are 
specified by initial velocity and initial spin orientation for deuterium and tritium. We apply 
the conservation of four momentum and conservation of intrinsic spin four-tensor. 
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II. CONSERVATION OF FOUR MOMENTUM 



The total four momentum before the coUision is equal to the total four momentum after 
the collision 

+Pf =Pr^+pf^ (2) 

where the superscripts d, t, n, and h label the four momentum vectors for deuterium, 
tritium, neutron, and helium, respectively. The subscript k — 1, 2, 3, 4 labels the space-time 
components of the four- vector. We follow the convention used by Synge and take the fourth 
component to be imaginary, thereby omitting the need to explicitly introduce a metric tensor 
for space-time. So the momentum four- vector has the form 

Pfe = (Pa,P4) = m(-^^a,^7) (3) 
c 

where 7 = (1 — {uaUa)/c^)~^^'^- Greek indices take values a — 1,2,3, and Latin indices 
always take values k — 1, 2, 3, 4. We use the convention that repeated indices are summed 
over their respective ranges. Here, Ua = dxa/dt are the ordinary three velocity components. 
Equation (2) is a conservation law that is valid in any inertial frame of reference. For our 
purposes, we assume that the laboratory frame is an inertial frame. The conservation of 
four momentum leads to four equations. 



III. ANGULAR MOMENTUM TENSOR 



In four dimensional space-time, the total angular momentum is an antisymmetric second 
rank tensor, Hmn = + H!^^\ with six independent components, for m^n = 1,2,3,4, 
where H!^^ is the orbital angular momentum and i^^n" the spin angular momentum [5]. 
The spin is represented by a four tensor because in four dimensions (space-time) there is no 
four- vector (of angular momentum) that is dual to the tensor if^^. 

As discussed above, we assume that in the laboratory frame of reference the impact 
parameter for deuterium on tritium is zero, so there is zero orbital angular momentum in 
the collision, i/"^'^^ = 0. In what follows, we drop the superscript and we use Hmn to denote 
the intrinsic spin four-tensor, H^^. 

For a particle that is not experiencing any forces, the intrinsic spin four-tensor is conserved 
along the particle's world hne [5]. When a coUision occurs at a space-time point, the sum of 
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the intrinsic spin four-tensors before the coUision must be equal to the sum of these tensors 
after the coUision, since the particles do not experience forces before or after the collision. 
We assume that the range of forces between particles is neghgible, so that we take the 
collision to occur at a point. The intrinsic spin four-tensor can be expressed in terms of the 
particle's four velocity, 

Xn^{-Ua,i^) (4) 

c 

and a unit spin four- vector, s„, according to [5] 

where Vt is the magnitude of the intrinsic spin, the four- velocity satisfies A^A^ = —1 and 
the unit spin four-vector satisfies, 

SnSn = 1. (6) 

The Levi-Civita symbol, ekimm is antisymmetric with respect to interchange of adjacent 
indices, it satisfies €1234 = +1 and has value -1 for odd permutations of its argument, so that 

( (S4A3 - S3A4) (S2A4 - S4A2) (S3A2 - S2A3) \ 

(S3A4 - S4A3) (S4A1 - S1A4) (S1A3 - S3A1) 

(S4A2 - S2A4) (S1A4 - S4A1) (S2A1 - S1A2) 

\ (S2A3 - S3A2) (S3A1 - S1A3) (S1A2 - S2A1) / 

The intrinsic spin four-tensor has the following frame- independent invariant, 

-Hkl Hkl = (8) 
Furthermore, the intrinsic spin four-tensor is orthogonal to the four velocity 

Hmn An = (9) 

by virtue of its construction in Eq. (5). 

The four components of the unit spin four- vector, Sn, are related to the four- velocity of 
the particle, A^, by the orthogonality condition [5] 

AnSn = (10) 

Equations (9) and (10) state that the physical spin vector hes in the three-dimensional 
hypersurface that is the three-dimensional physical space. Equations (5 ) through (10) are 
tensor equations and hence they are valid in any frame of reference. 



Hki = ifl 



(7) 



The fourth component of the unit spin four- vector, S4 can be expressed in terms of the 
other three components by use of the orthogonahty condition in Eq. (10), 

S4^ -{siUi + S2U2 + S3U3) (11) 

Using Eq. (11) to ehminate S4 from Eq. (6), we have a relation between the three components 
of spin, Sa, and the three velocity, up, which must be satisfied 

SaSa^l + -^{Sf^Uf^f (12) 

If we define three-component quantities by s = (si, S2, S3) and u = [ui, U2, M3), then we can 
write suggestively that the quantity s is normalized according to 

s2 = l + l(s.u)2 (13) 

Equation (12) or (13) shows that the three spatial components s of the unit spin four-vector, 
Sn, satisfy a normalization relation that depends on the particle's three velocity components, 
u. When the particle is at rest, u = 0, the three components satisfy s-s = 1, as expected. In 
the low velociy limit, we can interpret the three component quantity, S — fls, as the classical 
spin vector, where is the magnitude of the particle's spin. Of course, the transformation 
properties of the three components, S, depend on velocity and hence do not transform as a 
true three- vector. 



IV. INTRINSIC SPIN CONSERVATION 

We consider the collision, T (d , it) ^He where the initial conditions are such that both 
the deuterium and tritium spins are anti-aligned. The magnitude of total spin for deuterium 
is assumed to be = 1 and the magnitude of spin for the tritium nucleus is taken to be 
Qt — 1/2. The products of this inelastic collision are a ^ife nucleus, whose total spin is 
flh = 0, and a neutron, whose total spin is f2„ = 1/2. The conservation of intrinsic spin 
before and after the collision is given by 

Tjd I TTt _ Trh _|_ Trn (^A\ 

mn ' mn mn ' mn v / 

where = because the spin of He is zero. Each of the intrinsic spin four-tensors in 
Eq. (14) satisfies Eq. (9) by virtue of their construction according to Eq. (5). Each spin 
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tensor Hmn is a 4 x 4 antisymmetric matrix. Therefore, Eqs. (14) constitute six independent 
equations that relate the spin orientations to the velocity components, before and after the 
collision. Equations (14) are tensor equations that are valid in any inertial frame of reference. 

V. INITIAL CONDITIONS AND SOLUTIONS 

We have a total of ten conservation equations, where four of them are the conservation 
of four momentum given by Eq. (2) and six of them are the conservation of intrinsic spin 
given by Eq. (14). 

In solving the kinematic problem of the deuterium-tritium collision, the magnitude of spin 
for each particle is known, as described above. The orientation of the spins and velocities of 
deuterium and tritium are initial conditions that can be selected arbitrarily. We solve for the 
velocity and spin orientation components of the neutron, and the velocity components of the 
^He, in the laboratory frame of reference. For initial conditions, we take the tritium to be 
stationary at the origin of coordinates with velocity ul^ — (0, 0, 0), and we take the deuterium 
to be moving along the -|-x-axis, with velocity components specified by — {u, 0, 0), where 
u > 0. We take the deuterium spin direction to be along the -l-z-axis, specified by the spin 
orientation three- vector as = (0, 0, 1) and we take the tritium spin orientation three- vector 
to be = (0. 0,-1). The deuterium is moving such that s • u = 0, so = 1, and the initial 
condition = (0, 0, 1) satisfies Eq.(13). (As described above, the magnitude of the spin for 
every particle is fixed.) We solve for the following quantities: velocity and spin orientation 
components of the neutron, — (li", Uy, u^) and = (s", s^, s"), respectively, and velocity 
of ^He, u'^ = {u'^,Uy,u'^). Note that we use subscripts x,y,z and 1,2,3 interchangeably to 
denote components. 

We have a total of ten nonlinear equations, given by Eq. (2) and (14). There are a total 
of nine unknowns: three velocity components for the neutron, u^, three velocity components 
for the ^He, and three spin-orientation vector components for the neutron spin, s^. A key 
point is that, before solving the ten Eqs. (2) and (14), we must eliminate from Eqs. (14) the 
fourth component of the spin-orientation vector for each particle using Eq. (11). In this way, 
only the three spatial components, s^, of the spin-orientation four-vector sj^, appear in the 
intrinsic spin tensors in Eqs. (14). In the solutions of our ten equations, the spin-orientation 
vector for the neutron, s*^"') — (s^,s^,s^), and the neutron velocity, u^, must still satisfy 
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FIG. 1: The angle of neutron as a function of deuterium kinetic energy in the laboratory frame of 
reference for solution 1. 
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FIG. 2: The angle of helium as a function of deuterium kinetic energy in the laboratory frame of 
reference for solution 1. 
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solution 2 
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Deuterium energy [MeVJ 

FIG. 3: The angle of neutron as a function of deuterium kinetic energy in the laboratory frame of 
reference for solution 2. 

Eq. (12) or (13). We have checked that our solutions satisfy Eq. (12). 

Having ehminated the components, sf, s\, s2, and S4, the ten Eqs. (2) and (14) are solved 
by expanding them in a Taylor series to second order in the velocity components, u"^, and 
M, which we assume are small quantities compared to the speed of light, c. At least second 
order in the velocities must be kept to allow for the transmutation of mass, i.e., inelastic 
collision process. We find that there are two solutions to these ten equations, and each 
solution has a unique velocity for the neutron and for the ^He nucleus. The two solutions, 
which we call solution 1 and solution 2, are related by reflection in the x — y plane, which 
is perpendicular to the deuterium spin, = (0,0, 1), and contains the deuterium velocity 
vector, u'^ = (m,0,0). Note that the signs of the velocity components for the two solutions 
are as follows: 

solution 1 : < > 0, = 0, > 0; < 0, u'^ = 0, < 
solution 2 : < > 0, = 0, < 0; < 0, = 0, > 

Furthermore, we solve numerically the ten nonlinear (unexpanded) equations, given by 
Eq. (2) and (14), and compare these numerical solutions to the analytic solutions. The two 
solutions agree very well. 
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solution 2 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Deuterium energy [MeVJ 

FIG. 4: The angle of helium as a function of deuterium kinetic energy in the laboratory frame of 
reference for solution 2. 

The solutions show that the emitted neutron and hehum have velocity vectors that lie 
in the x-z plane. Figures 1 and 2 show the angle of the emitted neutron and helium, 
respectively, as a function of incident deuterium kinetic energy, for solution 1. Figures 3 and 
4 show that same quantities for solution 2, which is related to solution 1 by reflection, as 
described above. The numerical data used in our calculations is given in Table I. 

At the maximum cross section, where the deuteron kinetic energy is 107 keV, the angle 
of the neutron is On = arctan(n3/ni) = 82.85 degrees and the angle of the helium is 6h = 
arctan(M3/Mi) = 89.88 degrees. 

VI. DISCUSSION 

Within the context of special relativity, we have applied the conservation of four mo- 
mentum and the conservation of intrinsic spin to a collision of deuterium with a stationary 
tritium atom, in the laboratory frame of reference. We have assumed that the deuterium and 
tritium nuclei have spins that are anti-aligned and have zero impact parameter, correspond- 
ing to zero orbital angular momentum in the initial state. We have verified, as intuitively 
expected, that there are unique (classical) solutions for neutron (and helium) velocities, i.e., 
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TABLE I: Numerical data taken from Refs.[6] and [7] 



deuterium mass nid 2.01410177785 u 

tritium mass mt 3.01604927767 u 

neutron mass 1.00866491574 u 

^Hc mass m/, 4.00260325415 u 

speed of light c 299792458 m/s 

MeV per u - 931.49402823303 MeV/u 



there are unique directions for the emitted neutron and helium. Specifically, for a given 
energy of incident deuteron, there are two solutions that are related by reflection symmetry, 
see the above discussion. We have obtained the angular dependence of the emitted neutron 
and helium as a function of the incident deuterium kinetic energy. Our results may be 
compared with a quantum calculation [4] given by Kulsrud et al. 

Prom our classical relativistic analysis, it seems clear that we could in principle relax 
the constraint of zero impact parameter, so that we could impose some amount of orbital 
angular momentum for the deuterium and helium. Within a classical context (as opposed 
to a quantum calculation) the amount of orbital angular momentum imposed in the initial 
conditions is arbitrary and is a continuous quantity (not a discrete quantity, as in quantum 
calculations) related to the impact parameter. As described earlier, assuming non-zero an- 
gular momentum in the initial state would allow us to treat the case where the deuterium 
and tritium spins are aligned, which is presumably the more important case. As described 
earlier, in such as case, we will have orbital angular momentum appear in the products of the 
collision (neutron and helium). The continuum of possible angular momentum in the initial 
state would then lead to a continuum of angular momentum in the possible final (product) 
states. So the neutrons would not have a unqiue direction, but instead, they would have a 
distribution of possible velocities. The assumed initial probability distributions of impact 
parameters would then lead to a probability distributions for neutron velocities and angles. 
As described earlier, this is beheved to be the more important case. We leave this for possible 
future work. 
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